Asymmetric electron energy sharing in strong-field double ionization of helium 
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With the classical three-dimensional ensemble model, we have investigated the microscopic recol- 
lision dynamics in nonsequential double ionization of helium by 800 nm laser pulses at 2.0 PW/ cm 2 . 
We demonstrate that the asymmetric energy sharing between the two electrons at recollision plays a 
decisive role in forming the experimentally observed V-shaped structure in the correlated longitudi- 
nal electron momentum spectrum at the high laser intensity [Phys. Rev. Lett. 99, 263003 (2007)]. 
This asymmetric energy sharing recollision leaves footprints on the transverse electron momentum 
spectra, which provide a new insight into the attosecond three-body interactions. 
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Nonsequential double ionization (NSDI) of atom in 
strong laser field has drawn extensive researches in the 
recent years because it provides a particular clear man- 
ner to study the electron-electron correlation, which is 
responsible for the structure and the evolution of large 
parts of our macroscopic world |l|, |2J. The measure- 
ments of the recoil ion momentum distributions [3, |J| , the 
electron energy distributions J5j, [6| , the correlated two- 
electron momentum spectra 21 l8| , as well as numerous 
theoretical calculations [9l4l2| have provided convincing 
evidences that strong-field NSDI occurs in favor of the 
classical recollision model 



of Ar [15[ has not been well explored. 

Recently, the high resolution and high statistics exper- 
iments on double ionization (DI) of helium have made 
a great progress in unveiling the microscopic recollision 
dynamics in NSDI. The finger-like structure in the cor- 
related longitudinal (in the direction parallel to the laser 
polarization) momentum distribution from NSDI of he- 
lium by a 800 nm, 4.5 xlO 14 W/cm 2 laser pulses indi- 



131 ) . According to this model, 



the first electron that tunnels out of the atom picks up 
energy from the laser field, and is driven back to its par- 
ent ion when the field reverses its direction, and transfers 
part of its energy to dislodge a second electron. Though 
the recollision model describes the NSDI process in a 
clear way, the details of recollision remain obscure. For 
instance, at intensities below the recollision threshold, 
the underlying dynamics for the intensity-independent 
5U P (Up is the ponderomotive energy) cutoff in the two- 
electron energy spectra 14j-|l6[ and the dominant back- 



cates backscattering at the nucleus upon recollision [17 1. 
At a higher intensity, 1.5 xlO 15 W/cm 2 , Rudenko et 
al observed a pronounced V-like shape of the correlated 
two-electron momentum distribution 18], which is inter- 
preted as a consequence of Coulomb repulsion and typical 
(e,2e) kinematics. Theoretical studies have demonstrated 
that at the relatively low laser intensity, both the nuclear 
13 . |20J and the final-state electron 



Coulomb attraction 
repulsion 



20j,|2l( contribute to this novel structure. How- 



to-back emission of the correlated electrons from NSDI 
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ever, at the relatively high laser intensity, the roles of 
final-state electron repulsion and nuclear attraction for 
the V-like shape have not been examined. It is question- 
able whether the responsible microscopic dynamics for 
the V-like shape at this high intensity is similar to that 
at the relatively low intensity. 

In this Letter, with the classical three-dimensional 
(3D) ensemble model [lj, \2%, 



we examine the micro- 



scopic recollision dynamics in NSDI of helium by a high 
intensity (2.0xl0 15 W/cm 2 ) laser pulse. We find that 
the pronounced V-like shape of the correlated electron 
momentum in the direction parallel to the laser polariza- 
tion is a consequence of the asymmetric electron energy 
sharing in the recollision process, whereas neither the nu- 
clear attraction nor the final-state electron repulsion con- 
tributes to the V-like shape. This is different from that at 
relatively low intensity, where both the nuclear Coulomb 
attraction and final-state electron repulsion play signifi- 
cant roles in forming the finger-like shape. By separat- 
ing the recolliding electron from the bound electron, we 
find that the transverse (in the direction perpendicular to 
the laser polarization) momentum spectra for these two 
groups of electrons peak at different momenta. This dif- 
ference is ascribed to the Coulomb focusing in the trans- 
verse direction when the electron moves away from the 
core and can be understood as a footprint of the asym- 
metric electron energy sharing at recollision. 



The 3D classical ensemble model is introduced in [L; 
and widely recognized as an useful approach in study- 
ing high-field double ionization. In this classical model, 
the evolution of the two-electron system is governed by 
the Newton's classical equations of motion (atomic units 
are used throughout this Letter unless stated otherwise): 
^r = -V[Ke(n) + V ee (r 1 ,r 2 )} - E(t), where the sub- 
script i is the label of the two electrons, and E(t) is 
the electric field, which is linearly polarized along the 
x axis and has a trapezoidal pulse shape with four- 
cycle turn on, six cycles at full strength, and four-cycle 



turn off. The potentials are V ne (ri) — —2/y/rf + a and 



14e(fi,f2) = l/\/( r i — f2) 2 + b, representing the ion- 
electron and electron-electron interactions, respectively. 
The soft parameter a is set to 0.75 to avoid autoioniza- 
tion and b is set to 0.01 



12, 



221. To obtain the initial 



value, the ensemble is populated starting from a classi- 
cally allowed position for the helium ground-state energy 
of -2.9035 a.u. The available kinetic energy is distributed 
between the two electrons randomly in momentum space, 



and then the electrons are allowed to evolve a sufficient 
long time in the absence of the laser field to obtain stable 
position and momentum distributions [16j |. Note that in 
the classical model the first electrons are ionized above 
the suppressed barrier and no tunneling ionization oc- 
curs. 

Figures 1(a) and 1(b) display the correlated electron 
momentum distributions in the direction parallel to the 
laser polarization, where the laser intensities are 5.0 
xlO 14 W/cm 2 and 2.0 xlO 15 W/cm 2 , respectively. At 
5.0 xlO 14 W/cm 2 , the experimental observed finger-like 
structure is not reproduced (Fig. 1(a)). This is because 
of the large soft parameter employed in our calculation, 
which shields the nuclear potential seriously. Previous 
studies have illustrated that the finger-like structure is 
able to be reproduced when the realistic Coulomb poten- 
tial or a soft en p otential with a smaller screening param- 
eter is used [ljibol] . 

At the relatively high intensity, the overall V-like shape 
in the correlated momentum distribution is obvious. In 
contrast to the previous experimental result |18| . a clus- 
ter of distribution around zero momentum is clearly seen. 
Back analysis reveals that these events correspond to the 
trajectories where DIs occur at the turn-on stage of the 
laser pulse. For the soft potential employed in this Let- 
ter, the potential energy well for the second electron is 



— 2/V0.75 ~ —2.3 a.u., which is lower than that of realis- 
tic helium. In the classical description, the first electron 
can get ionized more easily at the expense of leaving the 
second electron near the bottom of the potential well 23 1 . 



Thus the first electron can be ionized very early at the 
turn-on stage of the pulse, leading to recollision occurs 
at the turn-on stage. This effect results in an overesti- 
mated contribution from the turn-on stage of the laser 
pulse to DI. In order to overcome this deficiency and fo- 
cus our study on the high intensity regime, we artificially 
exclude the events in which DI occurs at the turn-on stage 
of the laser pulse, as shown in Fig. 1(c). The correlated 
electron momentum distribution agrees excellently well 



with the experiment [lSj and the V-like shape is obvious 
though a soft parameter as large as a=0.75 is employed. 
We also performed further calculations by changing the 
soft parameter a after the first ionization [la . Il9j , and no 
noticeable change has been found in the V-like shape. It 
implies that the nuclear attraction does not contribute 
to the V-like shape, which is different from that at the 
relatively low laser intensity 13, |20| • 

It has been confirmed that at the relatively low inten- 
sity, the final-state electron repulsion plays an important 
role for the finger-like shape of the correlated electron 



momentum distribution 



20, 



21| . In order to examine 



the role of final-state electron repulsion in forming the 
V-shape at the high intensity, we have performed an ad- 
ditional calculation, in which the final-state electron in- 



teraction V r ee (ri,r 2 ) = 1/ ' \J{r\ — r 2 ) 2 + b is replaced by 



V ee (r 1 ,r 2 ) = ex]}[-Xr b ]/r b , where r b = y/(n - r 2 ) 2 +b 
and A = 5.0 |2oj]. As shown in Fig. 1(d), the V-like 
shape is still clearly seen, and no noticeable difference is 
found when compared to Fig. 1 (c) . Thus it confirms that 
the V-like shape is not a consequence of the final-state 
electron repulsion at this high intensity. 

The analysis above illustrates that neither the nu- 
clear attraction nor the final-state electron repulsion con- 
tributes to the V-like shape in the correlated longitudinal 
electron momentum at the high laser intensity. In order 
to explore the responsible dynamics for the V-like shape, 
we take further advantage of back analysis ll( . Tracing 
the classical DI trajectories allows us easily to determine 
the recollision time and the energy exchange during rec- 
ollision. Here, the recollision time is defined to be the 
instant of the closest approach after the first department 
of one electron from the core 12 1 . 



In Figs. 2(a) and 2(b), we have segregated the tra- 
jectories shown in Fig. 1(c) according to the energy dif- 
ference of the two electrons at time 0.02T after recolli- 
sion (T is the laser period. We have changed the time 
from 0.02T to 0.05T and the conclusions below do not 
change with the variation of this time.). Figures 2(a) 



and 2(b) display the correlated longitudinal electron mo- 
mentum distributions of the trajectories where the en- 
ergy difference is larger and less than 2.0 a.u., respec- 
tively. It is clearly shown that the events are clustered 
on the main diagonal when the two electrons achieve sim- 
ilar energies at recollision (Fig. 2(b)). In contrast, the 
correlated electron momentum distribution exhibits dis- 
tinct off-diagonal features when asymmetric energy shar- 
ing (AES) occurs (Fig. 2(a)). Based on these results, we 
can conclude that the AES at recollision is the decisive 
reason for the V-like shape in the longitudinal electron 
momentum correlation at the high laser intensity. 

In order to further understand the AES at this high 
laser intensity, in Figs. 2(c) and 2(d) we present the 
counts of DI trajectories versus laser phase at recollision. 
Figures 2(c) and 2(d) correspond to the trajectories from 
Figs. 2(a) and 2(b), respectively. It is found that in Fig. 
2(c), where AES occurs, recollisions cluster around the 
zero crossing of the laser field, while in Fig. 2(d), rec- 
ollisions occur close to the extremum of the field. Ac- 
cording to the simple- man model 13j , the electrons with 



the maximal recolliding energy return to the core near 
the zero crossing of the laser field. While these returning 
to the core near the extremum of the field possess lower 
recolliding energies. Figures 2(c) and 2(d) imply that 
the energetic recollisions often favor AES while the less 
energetic ones tend to have more symmetric energy shar- 
ing (SES). After distinguishing the recolliding electrons 
from the bound electrons we find that for 88% of the 
AES events (the events in Fig. 2(a)) the energy of the 
recolliding electron just after recollision is higher than 
that of the bound electron. It indicates that in the high 
rcturning-energy recollision, the recolliding electron only 
transfers a small part of its energy to the bound electron. 



24] , in which 



This issue is consistent with a recent study 
it has been demonstrated that the efficacy of electron- 
electron collisions decreases with the increasing collision 
energy. 

At the relatively low laser intensity, because of the 



lower recolliding energy, the efficacy of energy exchange 
at recollision is high. Thus AES is not serious and its 
contribution to the finger-like structure is negligible. At 
the high laser intensity, the low energy exchange efficacy 
at recollision makes AES play the dominant role in form- 
ing the V-like shape in the correlated electron momentum 
spectrum. Because of the dramatic AES, the two elec- 
trons leave the core with very different initial momenta 
and separate quickly. As a consequence, the final-state 
electron repulsion is weak and does not contribute to the 
V-like shape. 

More details of recollision can be obtained by inspect- 
ing the transverse momenta because the subtleties of the 
momentum exchange in the recollision process are not 
covered by the much larger momentum transfer taken 



from the laser field 



251 ]- In Figs. 3(a) and 3(b), we 



present the joint-probability distributions of the trans- 
verse momenta (along y axis) for the events shown in 
Figs. 2(a) and 2(b), respectively. Remarkably, in Fig. 
3(b) the distribution lies along the diagonal pi y +p2y = 0. 
This behavior indicates the strong repulsion in the trans- 
verse direction, which is in agreement with precious stud- 
ies 2a] . Contrarily, in Fig. 3(a) the population is clus- 
tered along the axes pi y = and pi y — 0, indicating 
different amplitudes of transverse momenta of the two 
electrons. This difference is more clear when separating 
the bound electrons from the recolliding ones. In the 
bottom of Fig. 3, we display the transverse momentum 
(Pjj_ = JPiy+Piz) spectra of the recolliding (red cir- 
cle) and the bound (black triangle) electrons separately, 
where Figs. 3(c) and 3(d) correspond to the events from 
Figs. 3(a) and 3(b), respectively. For the SES trajec- 
tories (Fig. 3(d)), the recolliding and the bound elec- 
trons exhibit similar transverse momentum distributions. 
Whereas for the AES ones (Fig. 3(c)), the difference 
in the distributions of the recolliding and bound elec- 
trons is remarkable: the spectrum of the bound electrons 
peaks near 0.2 a.u., while for the recolliding electrons the 
spectrum exhibits a maximum at 1.2 a.u. The different 



transverse momentum distributions for the SES and AES 
trajectories imply the different three-body interactions, 
which can be explored by monitoring the history of the 
DI events. 

We display two sample trajectories in Fig. 4. In the 
left column, the two electrons have equal energy after 
recollision (Fig. 4(a)), and achieve similar final longitu- 
dinal momentum (Fig. 4(c)). For the trajectory shown 
in the right column, the two electrons share unequal en- 
ergies upon recollision. The recolliding electron (solid 
red curve) obtains a higher energy at recollision (Fig. 
4(b)) but achieves a smaller final longitudinal momen- 
tum (Fig. 4(d)) due to the postcollision velocity }12j . 
The time evolution of the transverse momentum is more 
interesting. As shown in the bottom of Fig. 4, for both 
trajectories the two electrons obtain similar transverse 
momenta with opposite directions upon recollision. For 
the SES trajectories, both electrons experience a small 
sudden decrease in the transverse momenta just after rec- 
ollision (Fig. 4(e)). For the AES trajectory, the bound 
electron suffers a much larger sudden decrease in the 
transverse momentum while the transverse momentum 
of the recolliding electron does not change after recolli- 
sion (Fig. 4(f)). We ascribe the sudden decrease of the 
transverse momentum to the nuclear attraction in the 
transverse direction when the electron moves away from 
the core. For the SES trajectories, the two electrons leave 
the core with similar momentum, thus the nuclear at- 
traction plays a similar role in decreasing the transverse 
momentum, resulting in the distribution along the diag- 
onal piy + pi y — in Fig. 3(b). For the AES trajectory, 
the nucleus does not effect the transverse momentum of 
the recolliding electron because it leaves the core with a 
very fast initial momentum. While for the bound elec- 
tron, it takes a longer time to leave the effective area of 
the core due to the small initial momentum, leading to a 
significant decrease of the transverse momentum caused 
by nuclear attraction. The transverse momentum change 
of the electron is determined by Ap^ = J F±dt, where 



F± is the transverse force of the nuclear attraction. As- 
suming an electron that starts at a field zero near the 
region x=2 a.u. with initial momentum vj_ = 1.2 a.u. 
and evolves in the combined laser and Coulombic field, 
it takes a time about 10 a.u. for the nucleus to decrease 
v± to 0.2 a.u. 

Simply speaking, in the AES trajectory, because of the 
different initial momentum, the nuclear attraction plays 
different roles in "focusing" the transverse momenta of 
the bound and recolliding electrons, resulting in the mo- 
mentum distributions in Fig. 3(c). In other words, the 
different transverse momentum distributions of the recol- 
liding and bound electrons reflect the AES at recollision 
and provide a new insight into the attosecond three-body 
interactions. 

In conclusion, we have investigated the attosecond rec- 
ollision dynamics in NSDI of helium at 2.0 x 10 15 W/crn 2 . 
At the high intensity, the bound electron often shares a 
small part of the recolliding energy at recollision due to 
the low efficacy of energy exchange at the high recolliding 



energy. This asymmetric energy sharing is the decisive 
reason for the observed V-like shape in the correlated lon- 
gitudinal momentum spectrum at the high laser intensity. 
Because of the asymmetric energy sharing recollision, the 
bound electron leaves the core with a small initial mo- 
mentum. Thus its transverse momentum is strongly fo- 
cused by the nuclear attraction when it moves away from 
the core. Whereas the recolliding electron leaves the core 
so fast that its transverse momentum is not effected by 
the nuclear attraction. The different transverse momen- 
tum spectra of the recolliding and bound electrons act 
as a signature of the asymmetric energy sharing at rec- 
ollision and provide a new insight into the attosecond 
three-body dynamics. 
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FIG. 1: (color online) Correlated longitudinal electron mo- 
mentum distributions for NSDI of helium by 800 nm laser 
pulses. The intensities are (a) 0.5 PW/cm 2 and (b)-(d) 2.0 
PW/cm 2 . In (c) and (d), the trajectories where DI occurs at 
the turn-on stage of the trapezoidal pulse are excluded. In 
(d), the final-state e-e repulsion is neglected by replacing the 
soft Coulomb repulsion with Yukawa potential (see text for 
detail). The ensemble sizes are 2 millions. 
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FIG. 2: (color online) Correlated longitudinal electron mo- 
mentum distributions for the trajectories where the energy 
difference at time 0.02T after recollision is (a) larger than 2 
a.u. and (b) smaller than 2 a.u. (c)(d) Counts of DI trajecto- 
ries versus laser phase at recollision for the events in (a) and 
(b), respectively. The solid green curves represent laser fields. 
In all plots, the events where DI occurs at the turn-on stage 
of the trapezoidal pulse are excluded. 
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FIG. 3: (color online) (a)(b) Joint-probability distributions 
(on logarithmic scale) of the transverse momenta (along y 
axis) for the trajectories from Figs. 2(a) and 2(b), respec- 
tively, (c) Transverse momentum spectra of recolliding (red 
cycle) and bound (black triangle) electrons for the trajectories 
from (a), (d) The same as (c) but for the trajectories from 
(b). 
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FIG. 4: (color online) Two sample trajectories selected from 
Fig. 2(a) (right column) and Fig. 2(b) (left column), respec- 
tively. The upper, middle and bottom rows show the energy, 
longitudinal momentum and transverse momentum (along y 
axis) versus time for each electron, respectively. The energy 
exchange at recollision is clearly visible in the insets of (a) 
and (b). 



